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Here we propose to exploit the low energy bandwidth, small wavelength and penetration power 
of ultrashort pulses from XFELs for resonant Small Angle Scattering (SAXS) on plasma structures 
in laser excited plasmas. Small angle scattering allows to detect nanoscale density fluctuations in 
forward scattering direction. Typically, the SAXS signal from laser excited plasmas is expected 
to be dominated by the free electron distribution. We propose that the ionic scattering signal 
becomes visible when the X-ray energy is in resonance with an electron transition between two 
bound states (Resonant coherent X-ray diffraction, RCXD). In this case the scattering cross-section 
dramatically increases so that the signal of X-ray scattering from ions silhouettes against the free 
electron scattering background which allows to measure the opacity and derived quantities with 
high spatial and temporal resolution, being fundamentally limited only by the X-ray wavelength 
and timing. Deriving quantities such as ion spatial distribution, charge state distribution and 
plasma temperature with such high spatial and temporal resolution will make a vast number of 
processes in shortpulse laser-solid interaction accessible for direct experimental observation e.g. 
hole-boring and shock propagation, filamentation and instability dynamics, electron transport, 
heating and ultrafast ionization dynamics. 


PACS numbers: 82.53.Kp,52.25.Jm,52.38.Dx,52.38.Kd 


INTRODUCTION 


One of the essential elements in the generation of tran¬ 
sient hot solid-density plasmas by ultra-high intensit 


nsity 

ffl 


(UHI) lasers is the relativistic electron generation 
and transport dynamics 00. Near the laser focus, the 
current density of the relativistic electrons can exceed 
10 13 A/cm 2 j3|. This is intimately intertwined with the 
ionization dynamics and the complex evolution of the 
bulk return currents [g|, 9]. These are important due to 
the strong magnetic and electric fields generated at these 
current densities, up to 10 5 T and 10 14 V/m, as well as 
the rapid temporal and spatial evolution of the bulk tem¬ 
perature, ionization state, and hence resistivity by virtue 
of the electron-ion collision frequency, and anomalous re¬ 


sistivity from strong fields [5hl 
At present, a predictive understanding of high-intensity 
laser-matter interactions is severely hampered by the lack 
of self-consistent models for the ionization and recom¬ 
bination dynamics, coupled with the complex electron 
transport and collisions, and our inability to unravel 
this complexity with available experimental techniques 
in laser-only experiments. For example, effective ioniza¬ 
tion rates, scattering cross-sections and Ka and other 
self-radiation spectra are not precisely known in a non- 
thermal non-static dense plasma |jjJ. Additionally, nu¬ 
meric simulations are limited e.g. due to the short time- 
step necessary for kinetic simulation of fast plasma os¬ 


cillations, which can even reach classical electron orbit 
rotation times in an ion 0. Due to the large particle 
numbers in solid density plasmas the use of a full quan¬ 
tum mechanical ionization description will surely remain 
a dream for the next decades. 

With the availability of fourth-generation X-ray free 
electron laser (FEL) sources and their instrumenta¬ 
tion with ultra-intense short-pulse lasers many novel 
experiments can be envisioned that will e.g. enhance 
the diagnostic capabilities for femtosecond laser-solid 
interactions. Due to short duration of XFEL pulses 
and their penetration power through solids, coherent 
scattering allows probing of processes on femtosecond 
and sub-micron scale inside the solid plasma for the first 
time. Measuring the plasma opacity with high spatial 
and temporal resolution would allow for the first time 
to directly obtain experimentally the ion distribution 
during short-pulse ultra-intense laser interaction with 
solids and compare directly to simulations and ionization 
models used therein, or to study the plasma and hot 
electron dynamics through their influence on ionization 
states. 

We will show how the low energy bandwidth, small 
wavelength and penetration power of ultrashort pulses 
from XFELs can be used for Small-Angle Scattering 
(SAXS) on plasma structures 0 to not only scatter on 
free electrons but also to resonantly scatter at selected 
ions (resonant coherent X-ray diffraction, RCXD). 
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This allows to obtain information with high temporal 
resolution about the spatial structure of the distribution 
of electrons and ions (e.g. filament at ion, hole-boring, 
hot electron divergence and surface structures) as well 
as the energy distribution of free electrons (i.e. plasma 
temperature) [l34l6j with few nanometer resolution at 
the same time, as it is not possible with any other avail¬ 
able technique. It would allow unique tests of models 
e.g. for non-equilibrium ionization dynamics, heating, 
resistivity, electron transport or laser absorption. Pump- 
probe or pulse split-and-delay experiments would enable 
the dynamic measurement of the quantities, e.g. shock 
velocity, ion front and hole-boring velocity, or correlation 
spectroscopy. 

Resonant Elastic X-ray Diffraction or Anomalous 
Small-Angle X-ray Scattering UMi is a standard 
X-ray technique employing bound-free resonances 
that typically involve outer electron shells influenced 
by the surrounding structure of the atom to infer 
information on chemical or biological properties or 
physical structure (e.g of the lattice, spin distribution 
or magnetic domains) of a sample or molecule. In the 
present paper we discuss the possibilities that open up 
when transferring the methodology to bound-bound 
resonances in laser-generated plasmas in order to be 
sensitive to the ionic scattering signal and much less sus¬ 
ceptible to many-body phenomena and plasma screening. 


SCATTERING IN WARM PLASMAS 

The ionic scattering signal becomes visible against the 
non-resonant free electron scattering background when 
the X-ray energy is in resonance with an electron transi¬ 
tion between two bound states. In this case the scatter¬ 
ing cross-section dramatically increases due to the optical 
correction so that the signal of resonant X-ray scattering 
from ions silhouettes against the free electron scattering 
background despite the lower ion density. A schematic 
setup of a laser-plasma RCXD experiment is sketched in 
Fig. □ with the UHI laser and XFEL intersecting at the 
solid foil at approximately 90° to each other. Impor¬ 
tant parameters used in the following are summarized in 
Tab. U 

For a cold target the X-ray interaction is dominated 
by photoelectric absorption, see Fig. [2^ for the exam¬ 
ple of copper around the K-e dge. Once the temper¬ 
ature is sufficiently high and lower shells are not fully 
filled anymore due to excitation or ionization, channels 
for bound-bound transitions open up at photon energies 
just lower than the respective cold ionization edge, seen 
in Fig. [2^ for a temperature of 500 eV for Ka transi¬ 
tions in copper. When the XFEL energy Eq matches that 
of a bound-bound transition of an ion in a certain elec¬ 


target 



FIG. 1. Schematic setup of a RCXD experiment with a buried 
resonant layer target. Here the UHI is oriented normal to 
the target foil front surface and the XFEL penetrates the 
buried layer. The scattered photons are recorded far behind 
the target with a 2D areal detector after shielding and spectral 
filtering. 

UHI laser 

laser strength ag =6.8 
duration 7 - = 40 fs (FWHM) 

wavelength A = 0.8 pm 

wavelength A = 0.8 pm 

XFEL 

irradiated area A = 8 pm X 8 pm 
fluence lO^*-* photons/A 

^X F E L 8.290keV 

energy spread 20eV 

arrival time tarr = 8(3 fs past UHI 
duration typ. ion motion 

TABLE I. Summary of important parameters used in this 
paper. 


tronic configuration c, it will excite electron transitions 
between the two bound states accompanied by emission 
of X-rays. This can be treated as resonant scattering 
of the XFEL and be described by the complex valued 
change Ff es = F^(Eo) + iF'JjEo) of the ionic scattering 
form factor Fo = Fr = T [19| (inset in Fig. [2^). Here, 
r = Z — Q is the number of bound electrons of an ion 
with Z protons and total charge Q, the form factor fo of 
an electron is equal to unity since we consider only the 
limit of small scattering angles and we neglect the inter¬ 
nal electronic structure of the ions. 

For the more realistic case of finite XFEL FWHM band¬ 
width A E and a mixture of ions with different con¬ 
figurations F* es as found in warm plasmas, the aver¬ 
age optical corrections to the form factor of ions with 
T bound electrons are given by the sum over all con¬ 
figurations of a charge state weighted by their abun¬ 
dance n c (r ) and averaging over the XFEL energy profile, 
-Pf es = ^Ec inr (n c F* es (E)) E . The imaginary part of 
the optical correction caused by a single resonance can be 
derived from the product of its transition probability and 
the relative lower level population. The sum over all res¬ 
onances in all ions is proportional to the opacity r(Fo), 
t(E 0 ) = - 47rc g r 0 oM Er=i n r-Fr 0 (d is the foil depth 
and ro is the classical electron radius). Using Kramer- 
Kronig’s relation, the real part of the optical correction 
can be obtained. 


Plasma 


sim. box size 
cells/ pm 
material 
foil thickness 
preplasma 
particles per cell 10 ions 
(if fully ionized) 290 e- 


7.2 X 21.6 pm z 
250 

copper 

d = 1.6 pm in x-dir 
exp., scale 0.08 pm 
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Figure shows the opacity for the example of 
Copper for various plasma temperatures calculated by 
SCFLY 0 with distinct maxima visible. The typical 
naming scheme is following the resonant energy of ground 
state ions with T electrons, e.g. the He-, Li- and Be-like 
peaks for T = 2,3,4 respectively. Fig. [3] shows the con¬ 
tribution of specific ion charges to the total opacity for 
two plasma temperatures. For cold plasma, for example 
at the opacity Be-like peak around 8290 eV, ions with 4 
remaining electrons are responsible for less than 10 % of 
the opacity only. Due to the increase of abundance of 
highly charged ions at the expense of less charged ions, 
from 500 eV to 1000 eV temperature the contribution of 
ions with 5 electrons to the Be-like opacity peak decreases 
while that of ions with T = 4 and T = 3 increases. Yet, 
at no temperature does a single ion species dominate. 
In warm plasmas generally ions with different numbers 
of bound electrons contribute significantly to the total 
opacity at most XFEL energies. 

The reason is that spectator electrons in those shells not 
significantly screening the nuclear potential as seen by 
electrons participating in a transition do not influence 
its resonant energy, e.g. excited ions with specific T can 
have similar resonant energies as a ground state ion with 
electron number T — 1. To give a specific example, for Ka 
transitions an ion with T = 5 remaining electrons being 
in the excited state ls 2 2s 2 3p 1 with one spectator electron 
in the M-shell has generally similar resonant energies for 
transitions between the Is and 2 p !/ 2 and 2 p 3 / 2 levels as 
does an ion with T = 4 in the ground state ls 2 2s 2 . 
Therefore, in order to identify resonances, the number 
7 of electrons in the lower shells that do influence the 
transition energy of the respective resonance is a better 
parameter than T (in the above example the K and L 
shells). In analogy to F£ es we define 

£ (n c Frm E . a) 

Tl'y 

' c in 7 

extending the sum over all configurations with the same 
respective 7 regardless of the configuration of outer shells. 
Resonant scattering at a specific Eq can generally be at¬ 
tributed primarily to a single specific number of elec¬ 
trons in lower shells, justifying the labels of the opac¬ 
ity peaks given in Fig. [ 2)3 retrospectively with a slightly 
other meaning. E.g. from Fig. H]it can be seen that for 
the example of the Be-like peak around 8290 eV ions with 
7 = 4 dominate the resonant scattering largely indepen¬ 
dent of the plasma temperature. 


RESONANT SMALL ANGLE SCATTERING 
FROM ULTRASHORT LASER IRRADIATED 
SOLIDS 


The flux I(k x ,k y ) of scattered X-rays from scattering 
off electrons and ion resonances is given in first Born 




FIG. 2. Opacity of a copper plasma slab at different tem¬ 
peratures calculated by SCFLY & The inset in (a) shows 
exemplary the optical corrections for the specific configura¬ 
tion c = Is 2 2s 2 at its Ka 1 transition energy. The highlighted 
energy range indicates the XFEL energy and bandwidth used 
for the simulations in Fig. [5j (Note, the cold A-edge and 
Ao-resonances obtained from SCFLY exhibit a small shift 
compared to the NIST values [12], which however is not rele¬ 
vant for the qualitative morphology.) 


approximation by 


I(k X ,ky)= I(k) 

k z = 0 

H(r) = h iTee (r) + 


where!(k) oc H(k)H(k)* 
h hound (r) + h res (r ) 


( 2 ) 


= n free 


(r)fo 


z 

X 

r=i 


ni,r(r)Fr 


N-l 

X 

7=1 


n ii 7 (r)F 7 (r) 


where H(k) is the Fourier transform of A(r), assum¬ 
ing a quasi-static plasma during the XFEL pulse, XFEL 
pulse propagation along the z- axis and N is the max¬ 
imum occupation number up to the transition’s upper 
level. Near bound-bound resonances |A^ es (r)| increases 
rapidly and can reach values well exceeding 100 , then 
dominating over the non-resonant contributions. This 
is the reason for strong ion contrast with elemental and 
inner shell charge state specificity in the scattering sig¬ 
nal even at lower ion density than that of free electrons. 
Additionally the non-resonant free and bound electrons 
contribute dominantly to signal in other k regions than 
the resonant signal as the target keeps net neutrality on 
sufficiently large scales, and hence can generally be easily 
discriminated from the latter. 

A schematic setup of a laser-plasma RCXD experiment 
was shown in Fig. [T| with the UHI laser and XFEL in- 
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FIG. 3. Contribution of different ion species to the plasma 
opacity r oc nrF^' (solid lines) and nrF^ (dashed lines) for 
Copper ions with T = 3 to 7 bound electrons and plasma 
temperature 500 eV (a,b) and 1000 eV (c,d) from SCFLY. (b) 
and (d) show the relative magnitudes of the contribution to 
the opacity for different ion charges, averaged over 8280 eV — 
8300 eV. 



FIG. 4. Contribution of different ion configurations of the K 
and L shells to the plasma opacity r for Copper ions at plasma 
temperature 500 eV (a,b) and 1000 eV (c,d) from SCFLY cal¬ 
culations. Ions with 7 = 4 (red, orange) dominate regardless 
of plasma temperature. The highest specificity up to 1 keV 
plasma temperature is obtained at Eo = 8275 eV. Then, at 
500 eV, ions with 7 = 4 contribute almost 100%, at 1000 eV 
still ps 90% 


tersecting at the solid foil at approximately 90° to each 
other. In such a setup the XFEL penetrates the target 
foil in z-direction and is scattered at structures that are in 
the x — y— plane, e.g. filaments, shock fronts hole-boring 
compressions, plasma oscillations and others. The scat¬ 
tering pattern can be recorded behind the target with a 
2 D area photon detector. A simple experimental real¬ 


ization of the setup could be a flat foil with the XFEL 
beam intersecting the optical laser slightly less than nor¬ 
mal, penetrating the foil at a shallow angle. We note that 
in principle the presented RCXD method could also be 
used in a grazing incidence SAXS (GI-SAXS) geometry 
with X-ray penetration and scattering limited to the foil 
surface. The scattering pattern contains the full mode 
pattern of all electron correlations as well as that of res¬ 
onating ions. The real space distribution H(r) can even¬ 
tually be reconstructed using phase retrieval algorithms 
or b^employing more sophisticated holographic meth¬ 
ods j23l-[26j|. 

The scattering pattern 7(fe) does not only contain in¬ 
formation on the structure of spatial ion distribution but 
moreover information on distribution of charge states 
themself and plasma temperature. This information is 
encoded in the relative magnitude of real and imaginary 
part of the optical corrections, and it manifests itself in 
an asymmetry of the scattering pattern. If the XFEL 
energy Eo is tuned to match a resonance in ions with 
a specific 7 , also ions with more or less electrons in the 
lower shells may contribute as the resonances have a non¬ 
vanishing width (Fig. [4]). The contribution of such neigh¬ 
boring resonances manifests itself not only through the 
imaginary part of the optical correction, but also through 
the real part. Specifically, the real part of the optical 
correction at the center of a resonance vanishes, so any 
non-vanishing real part is due to neighboring resonances. 
Negative values of the optical correction real part are due 
to resonances at higher energy (smaller 7 ) while positive 
values are due to resonances at lower energy (larger 7 ) As 
the spatial distribution of ions with different values of 7 
is generally different, also the real part of the optical cor¬ 
rection will be distributed differently than the imaginary 
part, breaking the point-symmetry 7(fe) = /(—fe) of the 
scattering signal. This effect can be utilized by phase re¬ 
trieval methods to retrieve the imaginary and real parts 
separately and hence to infer the relative amount of ions 
with different charge state than that selected by Eq . This 
is especially interesting as different ionization models em¬ 
ployed in state-of-the-art simulation codes predict differ¬ 
ent distributions, thus allowing for a direct experimental 
validation. 

Furthermore, the ratio of real part to imaginary part 
of the optical correction, as well as relative amplitudes 
in a two-color XFEL experiment, can be used as a ther¬ 
mometer for the plasma temperature. The local electron 
temperature, or more generally the electron energy distri¬ 
bution, and its temporal evolution, determines the popu¬ 
lation of individual electron configurations and therefore 
the spectral shape of the optical corrections 14 
the number of electron holes in lower shells increases with 
plasma temperature, the real part of the optical correc¬ 
tion shifts from positive for small temperatures to nega¬ 
tive for higher temperatures. For a thermalized system 
this effect can be seen in Fig. 0] for the example of copper 
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ions, as simulated by SCFLY. The measurement of real 
and imaginary part of the optical correction can hence 
be used - in conjunction with an appropriate dynamic 
model - to determine the plasma temperature. 

So far, we discussed the change in the optical correc¬ 
tion real and imaginary part which has to be derived 
from the asymmetry induced into the scattering data via 
phase retrieval. However, using simultaneously multiple 
XFEL pulses with varying wavelength the exact spectral 
shape of the opacity curve and hence charge state 
distribution and temperature can be obtained directly 
experimentally. 


PIC simulations 


In the remainder of the paper we want to demonstrate 
the RCXD technique by example of simulation. We 
simulate RCXD images based on synthesized data from 
particle-in-cell (PIC) simulations of the laser interaction 
with a solid using the 2 dimensional collisional PIC code 
iPICLS 2 D [27| including field-ionization [H] and colli¬ 
sional ionization. We demonstrate the feasibility, sensi¬ 
tivity and physical interpretation of the method by com¬ 
parison of two different well established collisional ion¬ 
ization models, the Thomas-Fermi (TF) model (HJ and 
direct impact (DI) ionization 11]. In the TF model all 
ions are ionized up to the calculated local average charge 
state, thus producing a narrow distribution and neglect¬ 
ing the spread in charge states. The DI model is based 
on binary collisions generating a broad charge state dis¬ 
tribution but neglects recombination and only considers 
ground state ionization potentials; both models neglect 
internal excitations. We approximate the local distribu¬ 
tion of configurations (F^ and F") at the XFEL arrival 
time at the target foil t arr during post-processing with 
SCFLY, computing the temperature input-parameter us¬ 
ing the local average energy T e (r, t arr ) from the PIC sim¬ 
ulation with a time history proportional to that of the 
average over the simulation-box volume V, T e (r,t ) = 
Te(r, t arr ) (T e (r, t)) v / (T e (r, t arr )) v - The most impor¬ 
tant parameters used for the simulations are given in 
Tab. ID 

Since the PIC simulations are only 2 D in space we have 
to use a model to extend the plasma into the third dimen¬ 
sion (XFEL propagation direction z), in order to calcu¬ 
late the scattering signal. We define the target thickness 
in this direction to be D z = 1 pm but with only a cen¬ 
tral layer of d z = 0.2 pm thickness being resonant, while 
scattering in the remaining target volume only adds to 
the non-resonant background. Additionally we assume 
invariance of the plasma both in the resonant and non¬ 
resonant target volume in the direction of XFEL propa¬ 
gation. This setup corresponds to a sandwich target with 
the material of interest buried into material of neighbor¬ 


ing Z. The resonant layer and its thickness is chosen to 
be smaller than the typical structure size. The 2D sim¬ 
ulations are thus a good approximation of the real 3D 
situation while the non-resonant surrounding enables us 
to approximate the resonant vs. non-resonant signal ra¬ 
tios. 

Fig. [5] shows the simulated scattering patterns as they 
would be recorded by a detector with a pixel size of 15 pm 
at 1 m distance for the simulation employing the TF (a) 
and DI ionization model (c); Fig.[5fb) shows the scatter¬ 
ing pattern using the TF simulation but recalculating the 
ion distribution during post-processing by SCFLY based 
on the local plasma temperature and density in order to 
obtain more realistic and self-consistent ion distributions 
and ion configurations than TF alone. 

The lower left panels show the scattering signal taking 
into account only the free and bound electrons, i.e. no 
resonant transitions. This signal would be obtained far 
off any resonant energy and can clearly be separated in 
k -space from the resonant part. Since the electron distri¬ 
bution h free (r) + h bound (r) is rather smooth besides the 
sharp surface borders, the signal structure in transverse 
direction is weak and does not contain much useful infor¬ 
mation. 

The resonant signals (using only h re s in ©; lower right 
panels) and total signals show much more structure. This 
is, in this specific example, due to the filaments in ion dis¬ 
tribution and the ion front following the UHI laser pulse 
curvature, giving rise in the scattering pattern to trans¬ 
verse lines and a “bow-tie” like structure, respectively. 
Depending on the specific filament spatial frequency and 
width the scattering structure seen in transverse direc¬ 
tion changes (cp. Fig. [ 6 ]), which can be used to infer the 
electron dynamics that imprints on ionization } 30 ] as well 
as to benchmark ionization models. 

Additionally, a strong asymmetry is seen in the resonant 
scattering signal which is due to a difference in ion distri¬ 
bution for different charge states. As shown before, ions 
with differing 7 do not add significantly to the imagi¬ 
nary part of the optical correction. However, this is not 
strictly true for the real part. Hence, the real and imag¬ 
inary parts of H(k) are not only differing in amplitude 
(as they would if considering only one charge state and 
neglecting non-resonant electron scattering ), but also in 
phase. The strong asymmetry in the TF case is indicat¬ 
ing that ions with different 7 are spatially well separated 
while its absence in the case of DI is due to its broad 
charge distribution. 

By performing a phase retrieval on the scattering pat¬ 
tern, the real space distribution H(r) could be regained, 
individually for the real and imaginary part as shown 
in Fig |3 for the TF case (here directly taken from PIC 
simulation as used to generate the scattering patterns). 
Since the non-resonant signal corresponds mainly to that 
of a flat homogenous foil, its scattering signal is predom¬ 
inantly along the direction perpendicular to the foil sur- 
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FIG. 5. Simulated longitudinal charge state distribution along the laser axis (logarithmic scale) and 2 D average charge state 
from PIC simulations employing TF (a) and DI ionization (c) (target position is marked by dashed lines) together with the 
simulated scattering images (number of scattered photons per detector pixel) on a detector 1 m behind the target with pixel 
size 15 pm. Scattering signal for TF ionization with the ion distribution being constructed based on the plasma electron 
temperature by SCFLY during post-processing of PIC simulation is shown in (b). The signals as expected when the XFEL 
energy is tuned off resonances and the signals only considering resonant scattering are also shown for reference. Insets show 
the respective real space distribution of the resonant scattering centers h res (r), laser incident from left. The scattering wave 
vector scale bar in (a) applies to all scattering images. 



FIG. 6 . Scattering images for the PIC simulation employing TF ionization at four different XFEL energies corresponding to 
7 = 5 ( Exfel = 8240 eV), 7 = 6 (8202 eV), 7 = 7 (8162eV)and 7 = 8 (8124 eV). Insets show the respective real space 
distribution of the resonant scattering centers h res (r). All scales are the same as in Fig. [5] The scattering patterns illustrate 
from left to right the shift of ion distribution from compact distribution at the front to strongly filamented distribution following 
the electron filaments inside the foil. 


faces. All signal with a transverse component are hence 
due to resonant scattering, and as expected the optical 
correction’s real and imaginary part differ explaining the 
observed asymmetry. The real part exhibits the expected 
behavior in laser direction: at the foil surface the tem¬ 
perature is expected to be higher due to the laser irradi¬ 
ation, hence there should be more ions with inner shell 
holes (smaller 7 ) than inside the foil. With Eq centered 
at the resonance of ions with 7 = 4 this means more ions 
with 7 = 3 and hence negative real part of the optical 
correction at the surface and ions with 7 = 5 and positive 
real part inside the foil. Directly at the position of the 
ions with 7 = 4, the peak in the imaginary part, the real 
part is zero. 

To illustrate this we demonstrate how a simple poor- 
man’s phase retrieval can be used to derive the position, 
width and density of resonant ions and ions with one elec¬ 
tron more or less in the lower shells, as well as the local 
plasma temperature. We use the PIC data employing the 
TF ionization model at the resonant energy for 7 = 4, 
Eq = 8290 eV (Fig. [5^). We first approximate the real- 
and imaginary part of the scattering centers H(r) inside 
the target (cp. Fig. [7J) with box functions as shown in 
Fig. [H Computing the power spectrum for this model 
distribution we can fit the analytical result to the PIC 
synthesized scattering signal and obtain the density and 


position of resonant ions, as shown in Fig. [ 8 ] 

DISCUSSION AND OUTLOOK 

In this paper we have presented a novel diagnostic 
method that will enable for the first time to measure im¬ 
portant quantities in laser created plasmas overcoming 
limitations in spatial and temporal resolution of existing 
diagnostics. Tuning the XFEL energy to a specific energy 
selects a specific nuclear element with specific number of 
electrons in the inner atomic shells that will dominate the 
small angle scattering signal. Specificity depends on the 
XFEL energy and plasma temperature as the individual 
strength of resonances depends on the temperature de¬ 
pendent level populations but is generally possible to be 
better than 70% — 80%. The signals contain the full mode 
structure of ion-ion correlations that could be regained by 
numerical phase retrieval or holographic methods. The 
temperature sensitive variation of the opacity spectral 
shape can in turn be used to measure the spatial distri¬ 
bution of the instantaneous plasma temperature by using 
multiple-wavelength XFEL beams, and can hence be ap¬ 
plied for equation of state measurements. The presented 
method allows for unique measurements of the ion dis¬ 
tribution inside a warm or hot solid density plasma with 
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FIG. 7. Real space areal density of the target model multi¬ 
plied with form factor and optical corrections, h for the sim¬ 
ulation employing the TF ionization model (Fig. |5t]a)). (The 
absolute zero was shifted in each frame since only differences 
are relevant for the structure of the signal.) (a) h free + h bound , 
(b) 9ft(h res ), (c) 9ft(if) = h free + h bound + 9ft(h res ) and (d) 
$s(H) — 5s(h res ). (e-h) Respective lineouts along the laser 
beam axis. 
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FIG. 8. From the asymmetry of the scatter signal it is 
possible to infer the real and imaginary parts of H(r) via 
phase retrieval. For the simulation employing TF ionization 
(Fig. Oja)), already from a very simple model (red curve) ap¬ 
proximating 9ft (/i(z)) and 5s(h(z)) by rectangular shapes the 
width a, product of form factor and density and depth b, and 
position c of the resonant ion front can be obtained by fitting 
the power spectrum (eqn. under the right hand figure) to the 
simulated scattering data integrated over k y (ommitting the 
center three lines around k y — 0 containing the scattering of 
the bulk foil). 


high spatial and temporal resolution at the same time 
and hence can be used to directly test and improve ioniza¬ 
tion models on sub-micron few femtosecond scales for the 
first time. With refined simulations and theories - includ¬ 
ing e.g. self-consistently the ionization and excitation 
by laser accelerated electrons, plasma expansion, tempo¬ 
ral evolution of laser heating and cooling, recombination 


and de-excitation - one can derive the distribution of 
ion electronic configurations and hence construct the lo¬ 
cal history of ionization, electron energy distribution and 
density and eventually determine effects of non-thermal 
and transient electron dynamics on ionization dynamics. 
The spatial and temporal resolution are only limited by 
the XFEL wavelength and pulse duration. The temporal 
evolution of the properties may even be studied on fem¬ 
tosecond time scales employing more advanced schemes 
such as photon correlation spectroscopy using split and 
delayed pulses HIE). 

The presented method allows for the first time space 
and time resolved dynamical (via pump-probe or corre¬ 
lation spectroscopy) measurement on a few-nano meter, 
few-femtosecond level of many important properties of 
high-power laser generated plasmas such as the hole¬ 
boring, shock front and ionization front position, width 
and speed, mode structure of ion filaments and their re¬ 
spective growth rates and plasma temperature, poten¬ 
tially opening up a whole new field for understanding 
transient laser-plasma dynamics on a level not accessible 
to experiments before. 
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